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Abstract. We consider the continuity equation for open chaotic quantum 
systems in the semiclassical limit. First we explicitly calculate a semiclassical 
expansion for the probability current density using an expression based on 
classical trajectories. The current density is related to the survival probability 
via the continuity equation, and we show that this relation is satisfied within the 
semiclassical approximation to all orders. For this we develop recursion relation 
arguments which connect the trajectory structures involved for the survival 
probability, which travel from one point in the bulk to another, to those structures 
involved for the current density, which travel from the bulk to the lead. The 
current density can also be linked, via another continuity equation, to a correlation 
function of the scattering matrix whose semiclassical approximation is expressed 
in terms of trajectories that start and end in the lead. We also show that this 
continuity equation holds to all orders. 



PACS numbers: 03.65.Sq, 05.45.Mt 



1. Introduction 



The continuity equation is one of the basic concepts in physics, which appears 
in different contexts describing the conservation of a quantity. For example 
in electromagnetism it expresses charge conservation, in fluid dynamics mass 
conservation, and in quantum mechanics it represents the conservation of probability. 
Quantum mechanically the continuity equation can be expressed as 



where p(r,t) = \ijj(r,t)\ 2 is the probability density, 



j(r,i) = A[^( r , t )v^(r,t)-^(r,*)V^*(r,t)], (2) 



is the probability current density, and ip(r,t) is the solution of the time dependent 
Schrodinger equation with initial condition ip (r) at t = 0. A classical analogue of this 
continuity equation is [1] 




(3) 
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where 

<*.'.'>-«(-^). 

is known as the Morette- Van-Hove determinant, p = V r 5(r,r',t) is the momentum 
at r, and S(r, r',t) = J* pdq — Et is the action principal function, satisfying 
the Hamilton- Jacobi equation [2]. A first connection between the two continuity 
equations can be made through the WKB ansatz [3], which corresponds to taking 
^> WKB (r,i) = A(r, i)e* s ( r '*). By substituting this function into the Schrodinger 
equation, and neglecting higher order terms in H one arrives at the Hamilton-Jacobi 
equation, while the next order term in h leads to the classical continuity equation (3), 
upon the identification \ A\ 2 = \C\. 

It is then straight-forward, from this ansatz in the WKB approximation, to 
preserve the conservation of probability at the quantum level for small h, since this just 
relies on the classical preservation of probability. Moreover, this result is independent 
of the dynamics, and completely general for all Hamiltonian systems, though to our 
knowledge, there has not been much discussion about higher /l-dependent terms. 
The semiclassical methods which we use in this paper have the advantage though 
of providing us with additional information, allowing us for example to explicitly 
calculate quantum corrections to the survival probability. However, when applying 
these methods and semiclassical approximations, it is in fact an important issue to 
ensure that the unitarity of the quantum evolution is preserved, which in turn assures 
the fulfilment of (1). Our aim in this paper is to show explicitly that this is indeed 
the case for chaotic systems. 

In the context of transport through mesoscopic chaotic cavities [4, 5] this problem 
has been of wide interest. A possible way of describing the conductance of such a 
system is within the scattering approach [6]. The conductance is given in terms of 
asymptotic states in the leads projected onto Greens functions, which can be expressed 
semiclassically in terms of trajectories travelling from one lead to the other. The 
scattering matrix that describes the probability amplitudes of going from one channel 
to another, is unitary, which implies that the sum of all the probabilities should be 
equal to one. It has been possible to recover this result semiclassically by considering 
interference terms due to correlated trajectories, which were first considered for 
periodic orbit correlations in the context of spectral statistics. Moving beyond the 
diagonal approximation [7], the contribution of the first such correlated pair was 
found in [8], and this was later extended to all orders in [9, 10]. These ideas were 
then applied to transport in calculating the conductance, with the first contribution 
calculated in [11] and the extension to all orders in [12]. For the scattering matrix, it 
follows that the average sum of the probabilities is indeed one, and in fact there are 
also no fluctuations around this value [13, 14]. 

The continuity equation implies that the unitarity of the scattering matrix should 
be independent of the position of the cross-sections from which the scattering matrix 
is defined. Through this implication it is possible to show quantum mechanically 
that the scattering approach is equivalent [15, 16] to the Kubo-linear response theory 
[6], where the conductance can be written in terms of states inside the scattering 
region. It remains to be shown that this equivalence holds after applying semiclassical 
techniques. However, a mechanism for relating trajectories in the bulk with escaping 
trajectories was discussed in [17] for the continuity equation itself. They studied the 
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spatially integrated form of the continuity equation (1), which is given by 
d f 

+ JJ{r,t)-h x dx = 0, (5) 

where S is the cross-section of the opening and h x is the vector normal to this section 
at the point x in S. The survival probability p(t) is given by the integral J A dr p(r, t) 
over the volume A of the corresponding closed system. To satisfy this continuity 
equation it was shown [17] that it is necessary to include 'one- leg- loops', where the 
self-encounter now overlaps with the start or the end of the trajectory, among the 
correlated trajectories considered. Also, and further in [18], it was shown that 'one- 
leg-loops' were a basic ingredient needed to correctly calculate the survival probability. 
They are therefore necessary to recover unitarity, both of the survival probability and 
of the flow as expressed through the continuity equation. 

In this article we wish to build on these and previous semiclassical works to show 
that the conservation of probability described by the continuity equation is preserved, 
to all orders, when we take the semiclassical approximation to the propagator and 
perform the semiclassical expansion of the survival probability and current density. 
This approximation leads to expressions involving pairs of trajectories, and following 
the general philosophy of previous works we need to consider correlated pairs that 
can contribute in the semiclassical limit. As shown in [8, 9, 11, 12] these come from 
trajectories that have close self-encounters. For chaotic systems, the property of local 
hypcrbolicity means that we can construct a partner that differs in the encounter, 
leading to a small action difference and a contribution in the semiclassical limit, while 
the property of global ergodicity (and mixing) means that we can estimate the number 
of such pairs and calculate the contributions. These two properties, and such pairs of 
trajectories, are therefore responsible for the universal behaviour exhibited by chaotic 
systems, and the reason why we consider such systems here. 

In particular, among the correlated trajectories which we consider in this article, 
we also include those contributions coming from pairs of trajectories involving 'one- 
leg-loops' [17, 18]. Here we extend that work and show that the continuity equation 
is satisfied to all orders in the semiclassical approximation. For this purpose, we 
first re-examine the types of correlated trajectories which contribute and calculate a 
semiclassical expansion for the spatially integrated current density. By integrating 
the current density with respect to time we can compare with the previous result for 
the survival probability [18], to find agreement in line with the continuity equation. 
We then derive recursion relations for the possible trajectory structures that allow 
us to prove two things. The first is that, if we close the system, the semiclassical 
approximation preserves normalization as all higher order contributions cancel and 
the survival probability remains constant at 1. The second is a proof that the full 
continuity equation is satisfied to all orders. We also show that we can move via 
another continuity equation from the current density to a transport picture in terms 
of trajectories that start and end in the leads. All of these results are valid for time 
scales up to the Heisenberg time, and the direct extension to times longer than the 
Heisenberg time remains an open problem. We also restrict ourselves to a regime where 
the Ehrenfest time is much shorter than the average dwell time so that Ehrenfest time 
effects can be ignored. 

This article is organised as follows. We first introduce the semiclassical current 
density in section 2, where we also calculate the diagonal approximation, while in 
section 3 we extend the semiclassical calculation to higher orders. To begin our 
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proof that the continuity equation holds to all orders, we first shift to the energy 
domain via a Fourier transform in section 4 and recall the calculation of the survival 
probability. Our recursion relation arguments are presented in section 5, along with the 
proof that both unitarity and the continuity equation are satisfied in the semiclassical 
approximation. In section 6 we explore the connection to quantum transport and end 
with our conclusions in section 7. 



2. The semiclassical current density 

In order to calculate a semiclassical approximation to the current density, we start by 
writing the wavcfunction 



V(r,i) = J dr' if(r,r',t)^(r'), 



(6) 



in terms of the quantum propagator K(r, r', t) and the initial wavefunction ip { r ') at 
time t = 0. We then replace the exact quantum propagator with the semiclassical Van 
Vleck propagator [19] 

KSC ^ = J2^ £ ^(^>, (7) 

7(r'->r,t) 

where / is the dimension of the system, though in the following we will consider 
/ = 2. Sj(r, r',t) is the action along the path 7 connecting r' and r in time t, 

1 1 / 2 

and £> 7 = det — d 2 Sj(r, r', t)/drdr'\ cxp (— is the Van Vleck determinant 
including the phase due to the Morse index /i 7 . 

Substituting (6) and (7) into the expression for the current density (2), we obtain 



x 



I>7^e* (S, " S,,) tP7,f + P7',f] 



7(r'-»r,t) 
7'(r"^r,t) 

which involves a sum over pairs of trajectories 7 and 7', which have final momenta 
p 7) f and p 7 ' ; f at r. Due to the highly oscillating phase most of the contributions will 
cancel out upon averaging (for example over a local time average, though we do not 
include this explicitly in our notation) . The remaining systematic contributions come 
from pairs of trajectories that arc highly correlated and 'similar'. As in [17, 18], we 
neglect changes in the slowly varying prefactor and we expand the trajectories 7 and 
7', which go from r' and r" respectively to r in time t, around trajectories 7 and 7' 
which travel from the midpoint r Q = (r' + r")/2 to r also in a time t. We can then 
express (8) as 

r(M) = 4^ J dr ° £ ^ 7 <e^-^V(ro,P° r )P 77 <, (9) 

7>7'( r o— >r,t) 

where p° 7 , is the average initial momentum and p{ n , the average final momentum of 
the trajectories 7 and 7', and 

pw(r ,p) = |dq^ (r„ + fj (r Q - |) e"^P, (10) 
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is the Wigner transform of Vo( r o)> which arises from setting q = (r' — r"). This 
current density, or more accurately the integrated current density 

J(t) = [ j(r,t)-n x dx, (11) 
Js 

is the quantity we wish to evaluate semiclassically. The integral over the cross-section 
of the lead means that we are interested in trajectory pairs that start inside the system 
and end in the lead itself. 



2.1. Diagonal approximation 

The simplest semiclassical contribution to calculate is the diagonal approximation [7] 
where we pair the trajectories with themselves, 7 = 7'. Restricting ourselves to these 
pairs, (9) simplifies to 

j*ag (M) . ftx = _1_ |dr l^7lVw(ro,P 7 ,o)p f ,, 7 , (12) 

while the same treatment for the probability density leads to the similar result of 

/ iaS (M) = ^2 / dr ° E l^7lWr ,P 7 ,o). (13) 

7(r — >r,t) 

Performing the sum over trajectories in (12) using the open sum rule [11, 20, 21], 
the diagonal approximation becomes 

f^{v,t)-h x ^(^ C -^) , re S, (14) 



» W 1 r,p 

where w is the size of the opening, [i = l/r^ is the classical escape where Td is the 
dwell time, and (. . .) FiP is a phase space average. 

Integrating over the opening cross-section just leads to a factor w so that 

J dia s(t) = ^-1*, (15) 

where by supposing that the wave function has a well defined energy we can drop the 
average over phase space. By integrating with respect to time and setting p dlag (0) = 1, 
we obtain 

p dia s(f) = e - "*, (16) 

which is the classical decay for a chaotic system for long times. This result also follows 
directly from (13) by using the open sum rule and integrating over the volume of the 
system [17, 18]. 



3. Off-diagonal contributions 

To calculate higher order corrections we consider the contributions of trajectories that 
have close self-encounters. Highly correlated partner trajectories can then be found 
that follow the original trajectory almost exactly, but which differ in the encounter 
regions leading to a small action difference. An example of such a trajectory pair is 
given in Figure la. 

We are interested in calculating the contributions of trajectories with an arbitrary 
number of encounters, each of an arbitrary size, and we work along the same lines as 
[9, 10, 12, 13]. An encounter that involves I stretches of the trajectory is called an 
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Figure 1. A schematic picture of a trajectory 7 (solid line) that approaches itself 
twice in two self-encounters, and a partner trajectory 7' (dashed line) formed by 
crossing the encounters differently. In (a) the start of the trajectories is outside of 
the first encounter (211) while in (b) it has been shifted inside the encounter (111). 



/-encounter, so that the example in Figure la has two 2-encounters. Trajectories can 
be labelled by the vector v, whose elements vi list the number of /-encounters along 
the trajectory (and hence also along its partner). The total number of encounters 
is V = y^Vu while inbetween the encounters are long trajectory stretches called 
links. We make use of the fact that we can generate the possible configurations 
(or structures) of trajectories labelled by v from the related closed orbit structures 
(formed by connecting the start and the end of the trajectory). The number of links 
of the related closed orbit is L = ^ Ivi, and by cutting each of the links in turn (and 
moving the cut ends to the required positions) we generate the open trajectories. The 
number of trajectory structures N(v) with the same v is simply related to the number 
of closed periodic orbit structures (in fact for the closed orbits, we effectively count 
the same orbit L times and have to divide by that factor). When we cut a link of the 
closed periodic orbit, this generates two links in the trajectory structure, so that the 
total number of links is L + l. However, this usual contribution [12, 13], which are also 
called 'two-leg-loops' (211) in [17], is just one of the possible contributions. We can 
also shrink the link at the start of the trajectory so that the start moves into the first 
encounter, as in Figure lb. This case gives a different semiclassical contribution and 
corresponds to the 'one-leg-loops' (111) of [17]. We note that because the trajectory 
leaves the system at its end (which is placed in the lead) it cannot then return to a 
nearby point to have an encounter, and so we can only shrink the link at the start of 
the trajectory. Therefore the possibility that both the start and end point are inside 
encounters, the Oil case in [18], cannot happen. As such, the semiclassical contribution 
can be separated into these two cases: 

A where the start point is outside of the encounters (211), 
B where the start point is inside an encounter (111) 

3.1. Case A: two-leg-loops 

For structures corresponding to a vector v, the action difference between the 
trajectories 7 and 7' is given by AS — su in the linearized approximation. The vectors 
s and u contain the differences, along the stable and unstable manifold respectively, of 
the encounter stretches (of 7) in Poincare surfaces transverse to each encounter. Their 
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semiclassical contribution can be evaluated using an auxiliary weight, w v , A (s, u, t), of 
such encounters 

J„ )A (i) = fj,N(v) J dsduw v ^ A (s,u,t)c~ f * t ^ei su , (17) 

where N(v) is the number of trajectory structures corresponding to the vector v. The 
exponential e _ ^* exp is the average survival probability of the structures and involves a 
correction due to the proximity of encounter stretches during the encounters. Because 
of this proximity, if one encounter stretch survives then they all should do and so the 
exposure time t cxp involves only the time of a single traversal of each encounter. If 
we label the V encounters by a, with each involving l a encounter stretches that last 
t a nc , and we label the L + 1 links by i, with each lasting t if then the exposure time is 
simply 

L+1 v v 

icxp = ** + E Cc = t - E (*" - ^enc- ( 18 ) 
i— 1 a—1 a — 1 

Each encounter time is defined as the the time during which the stable and unstable 
separations of all the encounter stretches remain smaller than some classical constant 
c. They are therefore given by 

Cc = T ln 1 — r~ 1 — r> ( 19 ) 

A maxj \s a j x max j\u a j\ 

where A is the Lyapunov exponent, while s a and u a are vectors containing only the 
stable and unstable separations of the encounter a. 

The weight of encounters can be expressed in terms of an integral 

*.,A(.,«,t) - ^ n^vuA (20) 

where t enc is the total time the trajectory spends in the encounters, t cnc = ^^=1 'aCc 
The links should all have positive duration, and the weight is simply an L-fold integral 
over the first L link times, while the last link time is fixed by the total trajectory time 
t. When we perform the integrals the weight function becomes 

Co a. +\ ^ ~ ^» ^cnc) / 01 \ 

Wv As,u,t)= LlnL _ VUA . (21) 

To calculate the semiclassical contribution we will rewrite (17) as 

J vA (t) = nN(v) J ds du z v , A (s, m, t)e-^e^ su , (22) 

where z v ,a(s, u, t) is an augmented weight including the term from the survival 
probability correction of the encounters 

z vA (s, u, t) = w VtA {s, u, tJeE-C-- 1 )"*™ (23) 
„ (t - E a iM L U a (i + (la - l)MCc) 

Lw-vn a t" nc 

where we have expanded the exponent to first order. The semiclassical contribution 
comes from terms where the encounter times in the numerator cancel those in the 
denominator exactly [9, 10]. Keeping only those terms, we then obtain a factor of 
(2irh) L ^ v from the integrals over s and u. Finally we need the number of structures 
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Table 1. Contribution A of different types of trajectory pairs to the integrated 
current density, along with the number of structures for systems with time reversal 
symmetry. 

L V %^ NH_ 

(2)3 6 3 i^(-^ + ^-^ + iiL£) 41 
(2)1(3)! 5 2 H *Lgt (*3 - ^ + ^) 60 

(4) 1 4 1 £gl (-2*1 + 20 



AT(i?) corresponding to each vector, which are tabulated in [22], and we can obtain 
the result for trajectories described by the vector v of interest. 

As an example, in Table 1, we calculate the contribution of trajectory structures 
described by vectors v with L — V = 3, for systems with time reversal symmetry. We 
use a shorthand notation to describe the encounters in the vector v so that each term 
(l) Vl represents that the structure has vi /-encounters. These results can be multiplied 
by the number of structures and summed to give the third order correction for this 
case 

^^Vr + ^ ^r + ^r)- (24) 

3.2. Case B: one-leg-loops 

We can write this contribution as 

Jv,B(t) = V N ( V ) J ds du «»,b(s, u, i)e-^e^ SM . (25) 

Now that one encounter overlaps with the start of the trajectory, we have one link 
fewer [L in total) and an extra integral over the position of the encounter relative to the 
starting point. Starting with a closed periodic orbit (and dividing by the overcounting 
factor of L), we can cut each of the L links in turn and move the encounter on one 
side of the cut to the start. In total we obtain l a > copies of the same 111 involving 
the encounter a'. The augmented weight can then be expressed as a sum over the 
different possibilities, each of which involves an integral over the position of the start 
point inside the encounter, t a >, [17] 



ft ~ Sq ^enc) V 



z v , B (s,u,t) = £ l a , I " dt a , ^ ^TrA e^-D^c, (26) 

a , = 1 JO 11a W 

where the time of encounter a' is related to the starting position via 

Cc = ^' + yln j ,. (27) 

Because of the integrals over the position of the start point inside the encounter, the 
semiclassical contribution is calculated differently, using integrals of the type in [23] . 
With a change of variables however [18], the integral effectively gives a factor of t" nc , 
and we can write the augmented weight as 

, ra „ „ (Eg U° c ) ft Eg 'aO'' 1 Tig (1 + (ja ~ KJ ^ 

z v ,B(s,u,t) , (28) 



enc 
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Table 2. Contribution B of different types of trajectory pairs to the current 
density for systems with time reversal symmetry. 

L V %^ N(v) 

(2)3 6 3 J^(2g.-M£ + iW) 41 

(2)W 5 2 + 60 

(4) 1 4 1 H± ^(t) 20 



and treat it as before. 

Continuing with our example, in Table 2, we calculate this type of contribution 
for trajectories with L — V = 3, for systems with time reversal symmetry. These 
results can be multiplied by the number of structures and summed to give the third 
order correction for this case 

fie^ (20t 3 10/xt 4 41^ 



3.3. Results 



By simply adding the results for case A and case B, we can obtain the results for each 
vector and for each symmetry class. However, as can be seen in the example above, 
the Heisenberg time dependence involves only the value of L — V of the vector, and 
we can further sum over all vectors v with the same value of L — V to obtain that 
order in the expansion of J(t). The perturbative expansion is therefore in powers of 
t/Tn, where each term involves a finite expansion in powers of fit. We note that the 
same ordering of structures for transport quantities like the conductance leads to an 
expansion in powers of the inverse number of channels [12]. When we perform the 
expansion of the current density for the unitary case we obtain 

fie'^ ( fit 3 fi 2 ^ 



, , fie~^ ( fit 5 fi 2 t 6 7fiH 7 fiH 8 \ 



, / fit 7 401/i 2 * 8 643fiH 9 2fiH w 

Mt) - (-— + + (30c) 



ll^t 11 fiH 12 \ 



+ 



80640 322560, 
and zero for odd values of L — V . For the orthogonal case we obtain 



T2 l r 6 + 24 



J8(t) = TlT + l 24~ + ) ' (31C) 
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J4(t) = 2i — + — — — — + , (31d) 



^ 10 1 24 720 40320 , 

, , [ic-v* ( 16t 5 106/rf 6 121 M 2 t 7 61^ 3 i 8 , 01 . 

JsW = "j*" l - — + 15 30~ + -TO" (3l6) 



-H 

1321^ 4 t 9 2743m 5 * 10 
17280 + 1209600 
fic-^ (m 6 5419/zt 7 86801^ 2 t 8 53273/z 3 * 9 

V~3 420 + 10080 22680 ^ ^ 

1699^ 4 t 10 4063/i 5 * 11 55459/i 6 i 12 \ 



+ 5760 241920 + 159667200 



By direct comparison with the results obtained for the survival probability in [18] 
we can see that, term by term, 

§- t Pm(t) = -J m (t), (32) 

for both symmetry classes and all m < 6. Up to this order, it then follows trivially by 
summing over m that the continuity equation (5) is satisfied. We will prove that this 
result holds to all orders, and in order to do so we first consider the (inverse) Fourier 
transform of both the current density and the survival probability. 



4. Fourier transform 



To show that the continuity equation holds to all orders it is simpler to make a Fourier 
transform rather than showing, for each m, that all terms sum to zero for each power 
of t. In fact this is reminiscent of the situation for parametric correlations for systems 
without time reversal symmetry where it is simpler to show agreement with RMT (for 
times shorter than the Heiscnberg time) to all orders for the correlation function [24] 
rather than for the form factor [25]. 

For convenience we restrict ourselves to positive times, and we will consider the 
(one-sided) inverse Fourier transform of the current density J(t) 

POO 

J{w) = / dr J(rT H )e 2 ^, (33) 
Jo 

as well as of the survival probability p(t) 

POO 

P{uj)= / drp(TT H )e 2 ^ T , (34) 
Jo 

where r = t/Tn, and we will also use the number of open channels M = /uTh. In the 
Fourier space the continuity equation (5) becomes 

T H JH - {2ttilu)P(lu) = 1, (35) 

and this is the relation which we want to show semiclassically. We note that the 1 on 
the right hand side comes from the diagonal terms given in (15) and (16). 
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4-1. Transformed current density 

Interestingly the semiclassical contribution for the current density can be separated 
into a product of contributions over the encounters and the links of the trajectory 
in a similar way as for the conductance [12]. This means that we can obtain the 
semiclassical result very simply. For example for case A the contribution from 
trajectories with structures described by v is 



4 A H = -^ TT d^-W (36) 




where we have used the first expression for the exposure time in (18). We can 
perform these integrals as before, and because the Heisenberg times mutually cancel, 
we effectively get a factor of (M — 2tt\u;)~ 1 for each link and a factor of — (M — 2iriLdl a ) 
for each encounter. The contribution therefore simplifies to 

■^M = M*(«)(-l) y 1 (M_2^)^ • (3?) 

Considering case B, we have one link fewer (L in total) and one encounter, a' at 
the start of the trajectory pair. The integral over s a > and u a > is of course different 
from the integrals over the remaining encounters but we effectively obtain a factor 
of 1. We also remember that to obtain the structures for this case we started from 
the corresponding closed orbits, so that we divide by the overcounting factor i, and 
note that the encounter a' occurs at the start l a i times. Altogether, the contribution 
simplifies to 

J„,bH = - L ^ l « ' (38) 

which we can also write as 

j nN(v)(-l) v -i l a , \ njM - 2^l a ) 

J VtB {u) 1 ^ (m _ 2Wq/) J (M-2^ ' (39) 

to remove the restriction on the product over a. 
4-2. Transformed survival probability 

For the survival probability, as we saw in [18], because the trajectories that contribute 
start and end inside the cavity we additionally have the possibility that there is an 
encounter that overlaps with the end of the trajectory leading to three possible cases: 

A where the start and end points are outside of the encounters (211) , 

B where either the start or the end point is inside an encounter (111) and 

C where both the start and end point are inside encounters (Oil). 

We can again write the semiclassical contribution as a product of contributions 
from the links and encounters and for the first two cases we obtain very similar results 
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as for the integrated current density in equations (37) and (39). For case A we arrive 
at 

P.,aM = (40) 
p 

For case B, because the end of the trajectory is inside the system it can also be during 
an encounter. We therefore obtain an additional factor of 2 compared to the integrated 
current density and the result is then 

P.,bM = 2^M. (41) 
A* 

For case C things are slightly more complicated as we now have both ends of our 
trajectories inside different encounters and again one link fewer. The two encounters 
at the ends effectively give factors of 1 while the remaining encounters and links give 
their usual contributions. We also need to know the number of possible structures. 
Starting again with the corresponding closed periodic orbit structures described by 
the same vector v (and dividing by the overcounting factor L) we can generate the 
required trajectories by cutting the links which connect two different encounters a 1 
and (3'. We can then shrink the cut ends and shift them to the correct positions to 
obtain the trajectory structures. If we record in the matrix elements Af a ,p{v) the 
number of ways of cutting links that connect encounters a and (3 in all the periodic 
orbit structures described by v then we can write the contribution as 



V - CK ' \ j-^(M -2iriLul a ){M -2iriLul fj ) } (M — 2tt'iuj) l ~ 1 



5. Recursion relations 



In the limit p — ► (or — > oo) the quantum wavepacket remains inside the cavity 
and the survival probability is identically 1 for all times. This term comes from the 
diagonal approximation, which means that all higher order terms from correlated 
trajectories should vanish. As was shown for the survival probability, p(t), in [17, 18] 
one-leg-loops (cases B and C) are necessary to ensure this unitarity and cancel the 
terms from two-leg-loops (case A) which do not vanish when the system is closed. We 
first show that this holds to all orders by developing recursion relations between the 
different types of contribution. 



5.1. Unitarity 

By considering the sum over vectors v which share the same value of L — V = m, 
and summing over the different cases, we will show that P m (u}, p — 0) vanishes for all 
m > and both symmetry classes. This is equivalent to showing that the lowest order 
t term in the polynomial multiplying cxp(— pt) of the various p v (t) also sum to zero, 
so that p m (t) involves a polynomial with lowest term proportional to t m+1 . We start 
with the semiclassical results for each case when p = 0, and for case A we have 

P v , A (w, P = 0)= W7^ n v? + i a , (43) 
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while similarly for case B, the result is 
because J2 a 1 = V. For case C we obtain 

/ V (^ /i = ())= (V%_ I ■ (4 .-„ 




For this case, it is useful to rewrite the sum over a and (3 as a sum over the components 
of the vector v. M a ,fi{v) records the number of ways of cutting links, that connect 
encounters a and (3, in the periodic orbits described by v. However we can see that 
the important quantities are the sizes of encounters a and f3. Instead we record in 
Nk,i{v) the number of links that join an encounter of size fe to an encounter of size I. 

With this in place we can sum the different types of contributions, and then sum 
over vectors v with the same L — V = m 



(—l) mJ ~ 1 L-V=m 



p ^=°)=(&^t e (-D y n< 



V 



(46) 



Let us first consider the third contribution 



« a k,l 

To simplify this further we use the following recursion relation which can be deduced 
from [10, 22] . We have that J\fk,i (v) records the number of links that join a fc-encounter 
to an ^-encounter. Instead of cutting the link to make the trajectory structure, we 
imagine shrinking the link so that the k and Z-encounters merge to form a new (k+l— 1)- 
encountcr. By considering the number of ways that it is possible to shrink the link 
and create a smaller periodic orbit structure, we obtain the result 

A4» = (^-^-^W ^- 11 ), (48) 

where Vk+i-i is the (k+l — l)-th component of v and ^[W-^fc+J-i] \ s the vector obtained 
by decreasing the components Vk and vi by one and increasing the component Vk+i-i 
by one (so that Vk+i-i + 1 = v< t+h-i +l We also use the substitution 

N(v) = ( ~ 1)V _ n ° f ° JV(t;), (49) 
which allows us to rewrite (48) as 

JJ la ^«) = _ v i^+^N(v^ k + l -% (50) 

a 

so that (47) becomes 

p m ,o{^ = o) = -t^Si" E^E^rr- 11 ^^- 1 ]). (5i) 

V ' v k,l 
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As combining a k and ^-encounter reduces both L and V by one, the resulting vector 
v [k,l^k+l-i] g -(-jjj Yias the same value of L — V = m and by considering v [k,l^>k+l-i\ as 
a dummy variable v' it can be shown [10, 22] that 

L-V=m L-V=m 

E v^r- 1] N(v [k ' l ^ k+l - 1] ) = E -Ui-^iV). (52) 

v v' 

By identifying this dummy vector v' with v and substituting into (51), the total 
contribution to P m (w,^ = 0) in equation (46) simplifies to 



L-V=i 



(L-2V)-J2^k+i-i 



k,l 



N(v). (53) 



Concentrating on the sum over k and I, we define k' = k + 1 — 1 where fc' > ^ because 
k > 2. The sum then becomes 



E «fc+*-i = E E ^ = E - 2 H' = E ( fc ' - 2 H' 

fe,; ;>2 fc'>z fc'>3 fe'>2 

= ^ V - 2 ^ Vfc , = L - 2V, (54) 

fc'>2 fe'>2 

and so the term is square brackets in (53) is identically 0. 

This result shows that P m (oj, /i = 0) = and hence that p m (t, /j, = 0) = for all 
m > 0, for both symmetry classes. This is consistent with the fact that the survival 
probability should be identically 1 for a closed system (/i = 0). 

5.2. Continuity Equation 

Building on this, we are now able to treat the full continuity equation in the Fourier 
space. Again this will require re-expressing the contribution from the third case, and 
we also need to sum over different vectors v with the same value of L — V = m. If we 
define 

then the first two contributions to the transformed survival probability, see equations 
(40) and (41), can be written as 

L-V=m 

JVaM= E (Ml^) (*» 

L-V=m ~, 

P„„bM=- £ S is -^_»(„,M), (57) 

where we have replaced the sum over a by a sum over the components of the vector 
v. To simplify the third contribution from (42) we use a version of (50), modified for 
this situation to include the extra factors 

(-l) v MkAv) U a (M ~2nitol a ) 

(M — 2iriu)k)(M — 2iriujl) (M — 27riw) L_1 

(M - 27riw(fc + i - 1)) 1 ' J ' 1 ' 
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We can then rewrite the sum over the dummy vector v' = v [ k > l ^ k + l ~ 1 ] as a sum 0V er 
v and obtain 

L-V=ra 



E E 



(t + i-iK+i-i 

^ {M -2muj{k + l-l)) 



N(v,M). 



(59) 



We return to the continuity equation for the off-diagonal terms (J m (t) + 
dp m (t)/dt — 0) in the Fourier space (35). To ensure that the continuity equation 
holds, we have to check that 



M 



P m ,A{u) + 



(27Tiw) [P ro , A M + P ro ,B(w) + Pm,c(w)] = 0. 



(60) 



Writing the left hand side in terms of the sum over vectors, we have to evaluate 

L-V=m 



E 



lvi(M - Attiuj) 
(M - 2tt\ujI) 



+EE 



27ricjfc'wfe/ 
(M - 27riwfc') 



N{v,M), 



(61) 



i v 7 j fe'>; 

where fc' = fc + 1 — 1. Following similar reasoning to (54) we can simplify the double 
sum inside the square brackets as follows 



EE 



k'>i 



2-K\Lok'vy 

(M - 2TTiuk') 



E 



2irioj(l - 2)lvi 
(M - 2iriujl) : 



(62) 



so that (61) becomes 



L-V=i 

E 

V 



E 



lv ; (M - Attilo - 2niu){l - 2)) 



L-V=m 



- E 



E^ 



(M - 2-K\LOl) 

N(v,M) = 0, 



N(v,M) 



(63) 



since If; = L. This verifies equation (60) and shows that the semiclassical 
expansion satisfies the continuity equation for all m > 0. For the remaining diagonal 
terms (which can be thought of as corresponding to m = 0) this can be verified directly. 

Both of these proofs rely on being able to re-express the contribution from the 
third case in terms of a sum over vectors which has a similar form to the other 
two cases. This relation is then responsible for the fact that we can shift from the 
survival probability to the current density via the continuity equation and remove the 
possibility of having trajectories from case C (011). 



6. Implications for transport 

We have seen how by differentiating p(t), with respect to time, we obtain — J(t) in 
line with the continuity equation. For p(t) we have a picture involving trajectories 
that start and end inside the cavity, and we have three cases to consider. When we 
differentiate p(t), and shift to J(t), we arrive at a picture in terms of trajectories 
that start inside the cavity but end in the lead, effectively removing the third case 
(011) (and halving the contribution of 111). The next step is to repeat this process 
and differentiate again with respect to time. This leads to the more usual transport 
picture involving trajectories that start and end in the leads, where we can only have 
case A (211). 
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We can consider the conservation of the current density, which has its own 
continuity equation. For systems with constant potentials, like billiards, this 
continuity equation has the simple form <9j(r, t)/dt + V/(r,i) = 0, where f(r,t) is 
the second (antisymmetric spatial) derivative of the local density 

/(r, t) = ? -L ty* ( r , t ) V V(r, t) - 2 V^* (r, t) V^(r, t) + ^(r, i) V V (r, *)] . (64) 

To obtain a semiclassical approximation for this quantity, we express the wavefunction 
in terms of the semiclassical propagator using equations (6) and (7) and follow the 
same steps as in section 2 

/ SC M) = i^W / dr ° E ^ye*^-^'Wro,P° y )( P ; y ) 2 . (65) 

7,7'(r — >r,t) 

We are interested in the integrated version of this quantity F(t) = J drV 2 /(r, t), which 
can be expressed in terms of trajectories starting and ending at the lead cross-section. 
Instead of using (65), which involves trajectories which start and end inside the cavity, 
we first return to (64) and replace all derivatives with respect to r, originating from 
the derivative with respect to time of the current density, by derivatives with respect 
to r , which is defined before (9), and neglecting changes in the amplitudes. This 
replacement is possible in the case of energy-conserving dynamics. We can then rewrite 
the resulting expression as a divergence with respect to r and r Q in a similar way as 
it is done with respect to r to obtain the continuity equation itself. Upon applying 
Gauss' theorem for transforming the integrals with respect to r and r Q to surface 
integrals, and again supposing that the initial wave function has a well defined energy, 
we arrive at the semiclassical expression 

F (t) M / dxdx ' E ^>* (S7 ~ Sy W7'fV,77', (66) 

which is expressed as a sum over trajectories travelling from the lead to itself. We note 
that this result also gives the main semiclassical contribution for more general chaotic 
systems, and not only those with constant potentials. We can extend the connection 
to transport by projecting onto the channel basis 

n*)«^rE E D,D;,e^- s -<'\ (67) 

H a,b 7,7'(a— »6,t) 

where the sum over a and b is over the channels in the lead and then we sum over 
trajectories connecting these channels. The projection onto the channel basis changes 
the form of the stability amplitudes D 1 slightly, though they can still be treated using 
open sum rules [1 1] . This form of F(t) is the Fourier transform of a correlation function 
of scattering matrix elements. 

Because encounters cannot occur in the leads we only have a single case, 
corresponding to 211, and we can perform the sum over correlated trajectories using 
the open sum rule and an auxiliary weight function as before [12, 13]. However we 
have an additional contribution for systems with time reversal symmetry when the 
start and end channels coincide (a — b). Then we can also compare the trajectory 7 
with the time reversal of its partner 7' and we obtain a factor of 2 for this channel 
combination. This extra possibility corresponds to coherent backscattcring and must 
be considered more carefully when Ehrenfest time effects are important. The sum 
over channel combinations therefore gives a factor of M(M + k — 1), where k = 1 
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for systems without time reversal symmetry and k = 2 for systems with time reversal 
symmetry. The diagonal contribution is 

F» as W-M 2 (l + ^)e-^ (68) 
while we can simply express the contribution of trajectories described by a vector v 



as 



^J(t)+F(t) = 0, (70) 



^W=m(i + ^)^,aW. (69) 

We again shift to the Fourier space, where because of (69), the integrated 
continuity equation 
d_ 

dt' 

becomes, for the off-diagonal terms 

(M + k- 1) J A (a,) - (2niLu) [j A H + J B (w)] = - , (71) 

where the term on the right is what is leftover from the diagonal approximation for 
the orthogonal case. Rewriting the left hand side in terms of a sum over vectors we 
have to see if the following holds (dividing through by fi) 

lvi(2iriuj) (k — l)L 



E E 



E 



^ M ) = -(M^)- ™ 



(M - 2iriujl) (M - 2muj) 

Focusing on the left hand side, we recall that L = Ivi and so we can rewrite the 
first term in square brackets as 

_ ^ im(M-2wl) 
(M-2™I) ' 

so that we can combine it with the second term. We also separate the third term to 
rewrite the left hand side of (72) as 



I 

L-V= 



m v v 7 

For the unitary case, the second line vanishes, and because of the result (A. 6) in 
Appendix A, we can see that the sum in (74) is identically zero for each m and so 
the continuity equation is satisfied. For the orthogonal case, using the result (A. 9) 
in Appendix A, we can see that the terms in the sum in the first line for m = k 
cancel with the terms in the sum in the second line where m = k — 1. The only term 
remaining when we sum over all m is therefore the term from the first line where 
m=\. This corresponds to a vector with a single 2-encounter for which we can easily 
evaluate 

L ^ 1 ^lvi(M -2mw(Z-l)) ... 1 

E E \M-2I1) N ^ = - (M-2^y (75) 

as N(v) = 1. This term cancels exactly with the remaining term from the diagonal 
approximation, and we have verified (72) and hence (70) for both symmetry classes. 
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7. Conclusions 

The initial result of this article was the calculation of a semiclassical expansion of the 
integrated current density. The semiclassical approximation for this quantity involves 
correlated trajectory pairs which start at some point inside the system and end up 
escaping through the lead. The expansion was calculated by working along the lines 
of [13], but because the trajectories start inside the system we also need to include the 
contribution of 'one-leg-loop' diagrams introduced in [17]. These, and their extensions, 
were explored in detail in [18] for the survival probability itself, which semiclassically 
involves trajectories that both start and end in the bulk. Because of this, there are 
three cases that need to be considered for the survival probability rather than the 
two cases which exist for the current density, but both quantities can be connected 
through the continuity equation. This connection, however, does not trivially hold 
semiclassically, so we first showed directly that our calculation for the current density 
matches up with that for the survival probability in [18] which in turn was shown to 
agree with the supersymmetric random matrix results in [26, 27]. We note that our 
methods can also be applied directly to the quantities in the continuity equation (1), 
though we concentrated on the integrated version (5) here because of the interesting 
combinatorics. 

We then proceeded to show that, within the semiclassical approximation, this 
continuity equation is satisfied to all orders. This proof involved the recursion relation 
arguments presented in section 5 and is the main result of this article. This all hinges 
on our ability to re-express the contribution from the third case, where both the 
start and end point are inside encounters and which is only present for the survival 
probability, in terms of contributions related to the other two cases. The unitarity 
expressed by the continuity equation is therefore reflected in the combinatorial results 
derived from considering how many valid trajectory structures can be created by 
merging together encounters in more complicated trajectories. Moreover, we can 
expect this re-expression to lie behind any semiclassical situation where we move from 
a picture involving trajectories connecting points in the bulk to a picture in terms of 
trajectories connecting points in the lead. Interestingly, continuity can be thought of 
as a manifestation of the gauge (phase) invariance of quantum mechanics. Our proof 
of the continuity equation shows that this gauge symmetry would also be satisfied 
semiclassically, but the connection raises the intriguing question of whether treating 
the gauge invariance directly could lead to a simpler semiclassical proof and also shed 
more light on the combinatorial structure underlying this work. 

An important point, however, is the generality of the continuity equation. We 
have shown that it holds semiclassically for chaotic systems in the regime of times 
shorter than the Heisenberg time, and where Ehrcnfest time effects can be ignored. 
Moving directly beyond the Heisenberg time is probably the greater challenge, 
but Ehrcnfest time effects have provoked much interest recently, especially due to 
their appearance in real physical systems. They have been treated semiclassically 
[28, 29, 23], only for low orders (which are the most interesting physically), but this 
has yet to be generalised. One thing that we do know is that additional diagrams 
play a role with finite Ehrcnfest time, and these arc not included in the formalism 
used in this article (as they give no contribution in the limit we were working in). 
It is therefore an interesting challenge to generalise the extra lower order diagrams 
and completely extend the semiclassical treatment to the Ehrcnfest time regime. Of 
course, the continuity equation holds in all regimes and provides the perfect playing 
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ground to explore all the trajectory structures that contribute when Ehrenfest time 
effects are important, as well as requiring that 'one-leg-loops' and their generalisations 
are included. As well as Ehrenfest time effects, an obvious extension of the work 
presented here would be to include tunnel barriers at the lead. The possible back- 
reflection when trajectories try to escape leads to additional possible diagrams, as well 
as to other changes in the system [30] . 

The final result of this article was that the current density is itself related, via a 
continuity equation, to a transport quantity that is expressed in terms of trajectories 
that start and end in the lead. This result was also proved, in our semiclassical regime, 
using more recursion relations derived from [10, 22]. More importantly, it shows that 
we can, in a two-step process, move from a picture involving pairs of trajectories 
connecting points in the bulk to the more usual transport trajectories which connect 
points in the lead. From the survival probability, this involves differentiating twice 
(in line with the continuity equations) and taking a Fourier transform to arrive at 
a correlation function of the scattering matrix, which itself is closely linked to the 
conductance. 



Acknowledgments 

The authors wish to thank the DFG for funding under GRK 638 and FOR 760. 



Appendix A. More recursion relations 



To prove our second continuity equation we return to the combinatorics of the number 
of periodic orbit structures described by a vector v [10, 22], and we start with the 
unitary case. By considering the number of ways a 2-encounter could merge with 
a /c-encounter (to form a (k + l)-encounter) by shrinking the links connecting the 
2-encounters of the structure, they arrived at the relation 

*» N{V) = £ (fc + iX^+i + i) N{v i^ k+% {AA) 

k>2 

where ■y[ 2 < fc ^ fe + 1 ] i s the vector formed from v by combining a 2-encounter and a k- 
encountcr to form a (k + l)-encounter, so that Vi and Vk are reduced by one, while 
Vk+i is increased by one. We want to turn this relation into a version for open systems 
involving the extra factors we have, in a similar way as was done for parametric 
correlations in [24]. Including the extra terms and rearranging, we arrive at the 
following 

_ (k + l)v [ ^ k+1] (M - 2iriujk) . f9 , , , 

+ V- k+1 >-N(v [2 ' k ^ k+1] M) = 

r^ 2 (ilf-2*iu,(k+l)) " (V ' Mj 

where v^^ k+1 ^ = Vk+i + 1. Because this is identically zero, if we sum over all vectors 
with a common value of L — V = m the result is still zero 
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+ E '™ o_- n I Tvi -N(v [2 - k ^ k+1] , M) 



k>2 



(M - 2mu)(k + 1)) 



= 0. 



As combining a 2-encounter and a fc-encounter reduces both L and V by one, the 
resulting vector v l 2 - k ^ k + 1 ] has the same value of L — V = m as v. The important step 
is then to express the sum over the resulting vectors as 

L-V=m L-V=m 

J2 v^ k r k+ ^N(v^ k ^,M)= J2 v' k+ Mv',M). (A.4) 

V v' 

If we then identify the dummy vector v' with v in our original sum, we can rewrite 
(A.3) as 



L-V=rn 

E 



2v 2 (M - 2muj) + lvi(M - 2muj(l - 1)) 



(M - 4?riw) 



(>3 



(M - 2ttiujI) 



where I = k + 1. The first term can be included as the I 
so the result reduces to 

L-V=ra 



N(v,M) = 0, (A.5) 
2 term in the sum over I, 



E E 



l>2 



lvi(M - 2nkj(l - 1)) 
(M - 2iriujl) 



N(v,M) = 0. 



(A.6) 



For the orthogonal case, we can also create a valid periodic orbit structure if we 
shrink a link that connects a 2-encounter to itself so that the 2-encounter disappears. 
This means that there is an extra term in the recursion relation. Recasting the relation 
from [10, 22] into a form we require for our situation, we obtain 



2v 2 (M - 2tt\u) 
(M - Attiuj) 



(k + l)v [ ?*T k+1] (M - 2niu;k) „ . , 
N(v, M) + J2 ~ /J o_:../i. I Tvi -N(v^ k ^ k+1 ^ , M) 



+■ 



k>2 



(M - 2Triuj(k + 1)) 
iV>[ 2 ^l,M) =0, 



(A.7) 



(M - 2niLu) 

where v^ 2 ^ is the vector formed from v by removing a 2-encounter and L(v^- 2 ^) is 
the number of links that the new structure has. As well as the resummation in (A.4) 
we can also express 

L-V=m L-V=m-1 

L(v [2 ^ ] )N{v [2 ^\M) = LN(v',M), (A.8) 

V v' 

but now, because removing a 2-encounter reduces L by two and V by only one, the 
value of L — V of our new summation variable v' is one less than that of v. Using 
these resummations, and the fact that when we sum the relation (A.7) over all vectors 
v with the same value of L — V = m the sum is still zero, we obtain a result of 



L-V= 



E E 

v l>2 
L-V-m-1 



lv t (M - 2ttilu{1 - 1)) 



+ 



E 



(M - 2ttiuI) 
L 



N(v,M) 



(A.9) 



(M - 2?riw) 



N(v,M) = 0. 



This result, along with (A.6), allows us to prove (70), the second of our continuity 
equations for both symmetry classes. 



The semiclassical continuity equation for open chaotic systems 



21 



References 

[1] P. Gaspard 1998 Chaos, scattering and statistical mechanics Cambridge University Press, 
Cambridge 

[2] V. I. Arnold 1978 Mathematical methods of classical mechanics Springer, New York 

[3] M. Brack and R. K. Bhaduri 2003 Semiclassical physics Westview Press, Boulder 

[4] K. Richter 2000 Semiclassical theory of mesoscopic quantum systems Springer, Berlin 

[5] R. A. Jalabert 2000 In New directions in quantum chaos International School of Physics 

'Enrico Fermi' Edited by G. Casati, I. Guarneri and U. Smilansky IOS Press, Amsterdam 

143 145-222 

[6] S. Datta 1995 Electronic transport in mesoscopic systems Cambridge University Press, 
Cambridge 

[7] M. V. Berry 1985 Proc. Roy. Soc. A 400 229-251 

[8] M. Sieber and K. Richter 2001 Phys. Scr. T90 128-133 

[9] S. Miiller, S. Heusler, P. Braun, F. Haakc and A. Altland 2004 Phys. Rev. Lett. 93 014103 
[10] S. Miiller, S. Heusler, P. Braun, F. Haakc and A. Altland 2005 Phys. Rev. E 72 046207 
[11] K. Richter and M. Sieber 2002 Phys. Rev. Lett. 89 206801 

[12] S. Heusler, S. Miiller, P. Braun and F. Haake 2006 Phys. Rev. Lett. 96 066804 
[13] S. Miiller, S. Heusler, P. Braun and F. Haake 2007 New J. Phys. 9 12 
[14] J. Kuipers and M. Sieber 2008 Phys. Rev. E 77 046219 

[15] R. A. Jalabert, H. U. Baranger and A. D. Stone 1990 Phys. Rev. Lett. 65 2442-2445 

[16] H. U. Baranger, R. A. Jalabert and A. D. Stone 1993 Chaos 3 665-682 

[17] D. Waltner, M. Gutierrez, A. Gousscv and K. Richter 2008 Phys. Rev. Lett. 101 174101 

[18] M. Gutierrez, D. Waltner, J. Kuipers and K. Richter 2009 Phys. Rev. E 79 046212 

[19] M. C. Gutzwiller 1990 Chaos in classical and quantum mechanics Springer, New York 

[20] J. H. Hannay and A. M. Ozorio de Almeida 1984 J. Phys. A 17 3429-3440 

[21] M. Sieber 1999 J. Phys. A 32 7679-7689 

[22] S. Miiller 2005 Periodic- orbit approach to universality in quantum chaos PhD thesis, 

Universitat Duisburg-Essen 
[23] P. W. Brouwer and S. Rahav 2006 Phys. Rev. B 74 075322 

[24] T. Nagao, P. Braun, S. Miiller, K. Saito, S. Heusler and F. Haakc 2007 J. Phys. A 40 47-63 

[25] J. Kuipers and M. Sieber 2007 J. Phys. A 40 935-948 

[26] D. V. Savin and V. V. Sokolov 1997 Phys. Rev. E 56 R4911-R4913 

[27] D. V. Savin and H.-J. Sommers 2003 Phys. Rev. E 68 036211 

[28] S. Rahav and P. W. Brouwer 2006 Phys. Rev. Lett. 96 196804 

[29] Ph. Jacquod and R. S. Whitney 2006 Phys. Rev. B 73 195115 

[30] R. S. Whitney 2007 Phys. Rev. B 75 235404 



